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Supersymmetry breaking sectors generically have an approximate global U (1) _R-symmetry. Ther- 
mal effects tend to restore the _R-symmetry in the reheated early universe. We revisit the gravitino 
generation from the thermal plasma and we argue that an _R-symmetric phase suppresses the gaug- 
ino masses and thereby the production of the helicity ±1/2 gravitino component. For reasonable 
values of the hidden sector parameters the gravitino can account for the dark matter of the universe 
with a relic abundance characterized by a remarkable insensitivity to the reheating temperature. 
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Introduction. A generic feature of the supersymme- 
try breaking sectors is the presence of an exact or ap- 
proximate continuous global U{1) symmetry called R- 
symmetry. As pointed out at [1] a continuous R- 
symmetry is a necessary condition for supersymmetry 
breaking and a spontaneously broken i?-symmetry a suf- 
ficient condition if the superpotential is a generic func- 
tion of the fields. Hence, if there is a supersymmetry 
breaking vacuum an i?-symmetry, exact or approximate, 
must be generally expected [2]. Further support for the 
models characterized by an i?-symmetry comes from the 
requirement that the cosmological constant vanishes in 
the vacuum [3] or the necessity to explain large hierar- 
chies [4]. However, _R-symmetry must be broken in order 
that Majorana gaugino masses are generated. 

In the thermalized early universe the i?-symmetry, as 
we will discuss, although broken at low energies can be re- 
stored at high temperatures. This entails that the gaug- 
ino masses had been i?-suppressed; a fact that leads us 
to reconsider the relic gravitino abundance. Here, we will 
derive a new expression for the relic density of the ther- 
mally produced stable gravitinos which does not have the 
conventional linear dependence on the reheating temper- 
ature [5-9] . On the contrary, it depends on a temperature 
characteristic of the i?-symmetric (exact or approximate) 
thermal phase. This result implies a significant relax- 
ation of the reheating temperature bound in order that 
the dark matter constraint Vl^h? ^0.11 [10] is satisfied. 
Furthermore, it commends that the gravitino with mass 
in the range C'(MeV)-C'(GeV) can generically account for 
the dominant dark matter component of the universe. 

R-Symmetry and Gravitino Production. In gauge 
mediation [11] the gaugino masses m\ arise at one- loop 
from the i?- violating operator / d^OlnXW^Wa + h.c. 
where X the i?-charged spurion superfield that breaks 
supersymmetry. When the [/(l)ij symmetry is a symme- 
try of the vacuum soft masses for Majorana gauginos are 
not allowed; see [12] for a different scenario. 

A global U{l)fi symmetry broken at zero tempera- 
ture can, in principle, be restored at finite temperatures. 
The phase transition from an i?-violating phase to an R- 



preserving one has important implications on the thermal 
production rate for the helicity ±1/2 gravitinos: gauginos 
become massless hence, the reaction rate that generate 
gravitinos is not {mx/rriQ)^ enhanced but is simply given 
by the ^ T^/Mp factor. Let us call the temperature 
of the transition to the i?-symmetric phase Tq. Then 
for temperatures larger than Tq the production of grav- 
itinos is dominated by the helicity ±3/2 component i.e. 
for light gravitinos it is substantially eliminated. Only 
for temperatures T < Tq the helicity ±1/2 (Goldstino) 
production rate is recovered and, unless the ±3/2 com- 
ponent is overpopulated the yield Yq is dominated by the 
To and not from the reheating temperature Tj-h- 

The scalar soft masses that arise at two loops from 
the operator / d'^ein{X''X)Q''e^Q do not violate the R- 
symmetry. However, the contribution of sfermions to the 
Goldstino production is subleading whence, we can focus 
only at those 2 — 2 reactions that include a member of 
gauge supermultiplets and ignore those with chiral su- 
permultiplets altogether. 

Gravitinos can be also produced from the thermal 
plasma due to the top Yukawa coupling, an effect con- 
sidered in [9]. Like the gaugino masses, the A-terms re- 
quire interactions which violate the [/(l)/j symmetry and 
therefore, we expect them to be suppressed as well at high 
temperatures. 

Thermal Restoration of the U(1)r Symmetry. Let us 
consider the standard paradigm of gauge mediation: 



W = FX + XX(j)(j). 



(1) 



The (j) and arc the messenger fields. This superpo- 
tential yields universal soft masses, i.e. sfermions and 
gauginos have masses of the same order. The (1) is R- 
symmetric and a requirement for universal masses is the 
breaking of the i?-symmetry at the vacuum. At tree level 
the A"-direction is a pseudofiat direction with non-zero 
F-term. Everywhere along the pseudofiat direction the 
i?-symmetry is spontaneously broken apart from the ori- 
gin where the i?-symmetry is conserved. Pcrturbative 
corrections lift the flatness. The integrated out degrees 
of freedom also induce a correction SK to the canonical 
Kahler 
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6K = \X\yA^ 



(2) 
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where A the cut-ofF of the theory. Thereby, there is the 
quadratic term 4F^|Xp/A^ at the zero temperature po- 
tentiaL At finite temperature the messenger fields induce 
a thermal mass on the spurion of the form T'^\^\X\'^N / A 
where N are the complex degrees of freedom for the 
messenger fields (j) {<f>) considered to be in a fundamen- 
tal (anti-fundamental) representation of a GUT group. 
Hence, for models that exhibit spontaneous i?-symmctry 
breaking the critical temperature for the second order 
phase transition towards the broken R phase is [14] 



Tn = 



4 F 
N\A 



(3) 



Models with a Kahler potential of the form K = \X\'^ + 
\X\'^/A'^ — e\X\^/A^ where e a positive factor are some 
examples [13, 14]. This Kahler may originate from a mas- 
sive integrated out sector that breaks the i?-symmetry 
spontaneously. An example is the O'Raifeartaigh-type 
models proposed in [15] that contain fields with i?-charge 
i? 7^ 0, 2 or the model discussed in [16]. 

In a thermalized universe it is the messenger fields (/), (fi 
that control the thermal average value of the i?-charged 
spurion X-field and thus, the thermal restoration of the 
i?-symmetry. When the i?-symmetry is thermally re- 
stored the thermal average value is X{T) = 0. We can 
define the i?-breaking parameter bR{T) 



bR{T) 



Xn 



(4) 



For a second order phase transitions, like the one dis- 
cussed above, it takes the approximate values 6fl(r > 
To) = and bniT < To) = 1. 

Let us consider now theories that break the R- 
symmetry explicitly. For instance the superpotential (1) 
may be supplemented by a constant term dW = Wo that 
cancels the cosmological constant or an i?-violating mes- 
senger mass term 6W = M(j)(j). We can assume also here 
a corrected Kahler of the form K = \X\^ - \X\'^/A'^. The 
i?- violating terms result in a linear to X term at the zero 
temperature scalar potential for X that we write in the 
general form 

V ^Vo + ai{X + X^) + a2\X\'' + OiX^), (5) 

where a2 = 4F^/A^. The i?- violating vev of the X field is 
Xo = —axja-i- At finite temperature T < A the thermal 
average value for the i?-charged X field over the zero 
temperature vev Xo reads 



This is recast to 



X(T) 
Xo 



a2 + 4 A2r2 



bR{T) 



1 



(6) 



(7) 



where the Tq is given by the eq. (3). Here, at Tq there 
is no any second order phase transition. There is never a 



barrier and the symmetry breaking occurs via a graduate 
increase of the mean value of the X{T) field. At high tem- 
peratures the i?-symmetry tends to be restored i.e. hR{T) 
tends to zero. Indeed, for high enough temperatures the 
i?- violating terms Wo and M(p(j) are negligible and an ap- 
proximate i?-symmetry restoration takes place. We note 
that for the case of constant term Wo the zero tempera- 
ture vev is Xo = WqA^ /2FMp [17] and for the messen- 
ger mass case, after the translation X — s- M/\ — X, is 
Xo = M/X. Obviously, when T — )• the i?-symmetry 
breaking scale takes its maximum value, i.e. the zero 
temperature one, and when T — > oo the i?-symmetry is 
restored. In other words, the bfi(T) parameterizes the 
i?-symmetry breaking scale at finite temperature with 
respect to the zero temperature breaking scale. 

Gravitino Thermal Production Revisited. The ther- 
mal tendency to restore the i?-symmetry has important 
consequences. When the i?-symmetry is fully restored 
the gaugino masses vanish; for the case of approximate 
i?-symmetry restoration the gaugino masses are accord- 
ingly suppressed. The expected suppression is of the form 
m\{T) ~ bji{T)m\. Here we consider: 



mx{T) ^ X{T) 
mx Xo 



(8) 



where mx the zero temperature gaugino mass. Therefore, 
the gravitino cross section from the scattering processes 
off thermal radiations, (Tiqv), becomes thermally depen- 
dent. Considering the dominant QCD 2 — > 2 processes 
[7] the cross section is modulated by the bfj{T) factor 



^.(p,T)oc|^(^i + 6|(r)^^ 



(9) 



where 777^ the gluino mass and g the strong coupling. The 
gravitino yield is then given by 



dT 



{E^(T)v) rirad 

HT 



(10) 



see e.g. [6, 7]. The H is the Hubble parameter and 
nrad = C(3)r'^/7r2. Lightest Supersymmetric Particle 
(LSP) gravitino we can neglect the yield of the helicity 
±3/2 component for temperatures less than 



Trh < 10^^777^1 GcV^ , 



(11) 



Whence, focusing on the interactions of the helicity ±1/2 
modes the yield Yq is given by 



yg(r) = -D 



/v.(l/2) \ 



HT 




dTbliT) (12) 



where E^^^'' 



/(3?77|Af|,) and D the dilution fac- 
tor. We considered Y^{Trh) «C Yq{T) i.e. the dominant 
source of gravitino production are the scatterings in the 
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thermal plasma and any pre-inflationary abundance was 
diluted by inflation and 



5^(1/2) 



{T)v 



(T) 



(13) 



as we can see from (9). We note that the quantity in the 
brackets of (12) has no explicit temperature dependence. 

For models that exhibit exact thermal restoration of 
the f/(l)K symmetry the 6/j(T) can be approximated by 
a step function: 6fl(T) = for T > Tq and &fl(r) = 1 
for T < Tq. Hence, the gravitino yield (12) for Trh > To 
reads in this case: 



g*s{T) 

9*s{Tq) 



HT 



To (14) 



To 



where we took into account that T ^ Tq- The ratio 
g*siT)/ g^siTo) corresponds to the dilution factor where 
g*s the number of effectively massless degress of freedom 
[18]. Therefore, for T < 1 MeV, i.e. after nucleosyn- 
thesis, for a decoupled gravitino, fi ~ 100 GeV and for 
Trh > To the gravitino abundance is 



0.2 



To 



108 GeV 



GeV \ frngiii) 
J 



ITeV 



(15) 



i.e. it is controled by the Tq: the i?-critical temperature 
given by the eq. (3) . 

For models that break explicitly the i?-symmetry like 
those previously discussed the bii(T) is given by (7) and 
the relevant part of the integral (12) is 



dTbliT) 



2 



TqT 
T^ + T2 



Arctan 



(16) 

For reheating temperatures larger than the Tq, which 
are the cases that we are interested in and especially 
for Trh ^ To, the integral approximates to — To0,/i/2; 
where, 0rh = Arctan(rr?i/2o) a coefficient that takes val- 
ues: 7r/4 < 9rh < 7r/2. For reheating temperatures lower 
than To the integral (16) converges to {—)Trh value as 
expected. 

Therefore for Trh > To the (12) reads 



g*siT) 
g*siTo) 



^rad 



HT 



Y^^rh (17) 



To 



where we considered the coefficients (i.e. the g^s, g* fac- 
tors [18]) in front of the integral to be dominated by the 
value given at the temperature Tq. The resulting grav- 
itino abundance is 0rh/'2 times the (15) 



0.1 



rhTo 



108 GeV 



GeV 



ITeV / 



(18) 



The expression (18) has a very mild dependence on the 
reheating temperature. We can say that the gravitino 



abundance is reheating-tempcrature-indepcndent as long 
as the reheating temperature is larger than Tq and not 
so large that the 3/2 component to be overproduced. In- 
deed, for temperatures Trh > 10 Tq the 9rh takes values 
only between 1.50 and 1.57. Substituting 9rh ~ 1-5 and 
To = 4i^/(V^AA) where 4F ~ 
expression (18) attains the form 



16.6 X IO^^^^QjjV the 



n/,h^ ~ 0.15 X 



16.6 



IQio GeV^ 



AA 



ITeV 



(19) 



We mention, that in the case of approximate R- 
symmetric models the gravitino abundance is slightly 
smaller by a factor of 2/9rh compared with the case of 
exact restoration of the i?-symmetry. The expectation 
is that the gravitino production is stronger suppressed 
when an exact restoration of the i?-symmetry takes place 
than when the restoration is approximate. This small dis- 
crepancy originates from the fact that we simplified the 
thermal evolution of the X{T) by assuming a step-like 
behaviour while a second order phase transition is not 
a discontinuous process; instead at the critical temper- 
ature To the barrier vanishes and the transition occurs 
smoothly. The Tq value of the integral (14) and thus the 
Yq bound in that case is a conservative one. 

Discussion and Conclusions. The expressions (15) 
and (18), (19) are a remarkable result. Firstly, the grav- 
itino abundance exhibits an insensitivity to the grav- 
itino mass and the reheating temperature. Secondly, 
the quantities which control the yield are the supersym- 
mctry breaking fundamental parameters A and A. For 
\/N = 0{1 - 7) [11] and mg - 1 TeV, the gravitino does 
not overdose the universe if 



AA > 10^" GeV 



(20) 



When the lower bound is saturated it can account for 
the dominant dark matter component; for instance, when 
A ^ 10"^ and A = O(GUT) scale, the gravitino is the 
dark matter of the universe. It is interesting to note that 
a GUT scale cut-off is expected in many theories like 
those discussed here. This fact also implies the smallness 
of the coupling since in the IR theory the coupling A is 
expected to be suppressed, see e.g. [19]. 

In addition, considering that the supersymmetry 
breaking local minimum must be thermally preferred we 
conclude that small values for the coupling A < 10"'' 
are cosmologically favourable [20]. On the contrary, the 
increase of the coupling A value renders the supersym- 
metry preserving vacua more attractive. We thus find an 
interesting window of values where the supersymmetry 
breaking vacuum is thermally selected and the gravitino 
(nearly automatically) does not overdose the universe, 
or even accounts for the dominant dark matter compo- 
nent. This parameter space also specifies the gravitino 
mass range. 

The interactions of the hehcity ±3/2 that are T^/Mp 
suppressed impose an upper bound on the Trh- Accord- 
ing to the relation (11) and for AA > 10-^° GeV, a grav- 
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itino, for instance, with mass 1 GeV (100 MeV) con- 
strains the reheating temperature to be T^h ^ 10^^ GeV 
(10" GeV) which is about 10^ (10*^) times relaxed rela- 
tively to the conventional bound on the reheating tem- 
perature for the stable gravitino. For tUq < 1 MeV the 
yield of the 3/2 gravitino component can be safely ne- 
glected [21] and thus the total yield is given solely by 
the (14) or (17). An interesting remark is that when the 
conventional (i.e. ignoring the 6/{(T) factor) freeze-out 
gravitino temperature is larger that the Tq the light grav- 
itino never actually equilibrates. In other words, in this 
case there is no freeze-out temperaure for the gravitinos. 

The i?-symmetry is generally expected to be restored if 
large reheating temperatures were realized i.e. Trh > Tq. 
In the minimum of the finite temperature effective po- 
tential the entropy maximizes. Here, the thermalized 
degrees of freedom are expected to be the (M)SSM de- 
grees of freedom and the messenger fields. The later 
have tree level mass ~ XX and become light at the ori- 
gin. Thermalized hidden sector fields that become light 
at the X = Xq could render the _R-breaking minimum 
thermally favourable rather than the origin X ^ 0. How- 
ever, this requires i?-violating interactions. Furthermore, 
a thermalized hidden sectror that couples directly to the 
Goldstino would lead to a large driQ/dt with a Goldstino 
production rate scaling like T^/F^ [22]. If such a hid- 
den sector exists it must be nearly unpopulated. On the 
other hand, the messengers although thermalized give a 



suppressed Goldstino production rate that can be ne- 
glected since A ^ 1. These arguments indicate that an 
i?-symmetry restoration probably takes place in a high 
temperature environment of the early universe. 

The mx ^ F/ A is of the order of the electroweak scale 
thereby, it is the smallness of the Yukawa coupling at the 
messenger superpotential that makes the Tq large sug- 
gesting a 0(MeV)-C'(GcV) mass range gravitino for dark 
matter (18). 

The conclusion is that the gravitino cosmological prob- 
lem can be actually absent in the most general class of 
gauge mediation models without including additional in- 
gredients or assumptions. This is deduced from the ex- 
pected gaugino mass i?-suppression described by the eq. 
(8). Moreover, much like the neutralino LSP which has 
a fixed thermal abundance as long as the Trh is larger 
than the freeze-out temperature, so the gravitino LSP 
has a relic abundance fairly insensitive to high reheat- 
ing temperatures as long as T^h > Tq. Barring that this 
fact is good news for thermal leptogenesis we emphasize 
that the gravitino can be indeed considered a compelling 
candidate for the bulk dark matter of the universe. 
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